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Abstract 

The measurement of the nuclear structure function F2'{x, Q^) at the future electron-ion collider 
(EIC) will be of great relevance to understand the origin of the nuclear shadowing and to probe 
gluon saturation effects. Currently there are several phenomenological models, based on very dis- 
tinct approaches, which describe the scarce experimental data quite successfully. One of main 
uncertainties comes from the schemes used to include the effects associated to the multiple scatter- 
ings and to unitarize the cross section. In this paper we compare the predictions of three distinct 
unitarization schemes of the nuclear structure function which use the same theoretical input to 
describe the projectile- nucleon interaction. In particular, we consider as input the predictions of 
the Color Glass Condensate formalism, which reproduce the inclusive and diffractive ep HERA 
data. Our results demonstrate that the experimental analysis of will be able to discriminate 
between the unitarization schemes. 
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I. INTRODUCTION 



The measurement of the nuclear structure function in deep inelastic electron-nucleus 
scattering (DIS) is the best way to improve our knowledge of the nuclear parton distributions 
and QCD dynamics in the high energy regime (See, e.g.[l, 2|). However, after more than 30 
years of experimental and theoretical studies, a standard picture of nuclear modifications of 
structure functions and parton densities has not yet emerged. Fixed target DIS measurement 
on nuclei revealed that the ratio of nuclear to nucleon structure functions (normalized by 
the atomic mass number) is significantly different from unity. In particular, these data 
demonstrate an intricate behavior, with the ratio being less than one at large x (the EMC 
effect) and at small x (shadowing) and larger than one for x ~ 10~^ (antishadowing). The 
existing data were taken at lower energies and therefore the perturbative QCD regime 
(Q^ > 1 GeV^) was explored only for relatively large values of the (Bjorken) x variable 
{x > 10~^). Experimentally, this situation will hopefully change with a future high energy 
electron- ion collider (EIC) (For recent reviews see, e.g. [4, 0]), which is supposed to take 
data at higher energies and explore the region of small x {x < 10~^) in the perturbative 
QCD regime. 

The theory of nuclear effects in DIS is still far from being concluded. The straightforward 
use of nucleon parton distributions evolved with DGLAP equations and corrected with a 
nuclear modification factor determined by fitting the existing data as in Refs. 6l4ll| is well 
justified only in the large region and not too small x. Moreover, these approaches do 
not address the fundamental problem of the origin of the nuclear shadowing and cannot 
be extended to small x, where we expect to see new interesting physics related to the non- 
linear aspects of QCD and gluon saturation (For reviews see Ref. [l^). Currently, there 
are several phenomenological models which predict different magnitudes for the shadowing 
in the nuclear structure function based on distinct treatments for the multiple scatterings of 
the partonic component of the virtual photon, assumed in general to be a quark-antiquark 



(qq) color dipole. Some works 
the Glauber- Gribov formalism 



16] address the origin of the nuclear shadowing through 
171,1181 in the totally coherent limit (/c ~ l/2mj<qx ^ Ra, 



where Ic is coherence length), which considers the multiple scattering of the color dipole with 
a nucleus made of nucleons whose binding energy is neglected. In the high energy limit, the 
eikonal approximation is assumed, with the dipole keeping a fixed size during the scattering 
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process. In this approach the total photon-nucleus cross section is given by 



where \ip{r, z)\'^ is the probability of the photon to split into a qq pair of size r and adA{x, r) 



is the dipole-nucleus cross section, which is expressed as 
adAix,r) = J d%2 



13| 



1 -exp ( --]^ATA{h)adp{x,r] 



(2) 



with T^(6) being the nuclear thickness function and (Jdp(x,r) is the dipole-proton cross 
section. It must be stressed that once (Trfp(x, r) is fixed, the extension to the nuclear case 
is essentially parameter free in this approach. In the Glauber formula it is assumed 
that the dipole undergoes several elastic scatterings on the target. Although reasonable and 
phenomenologically successful this assumption deserves further investigation. This model 
can be derived in the classical approach of the Color Glass Condensate formalism [l^. 

Another approach largely used in the literature is based on the connection between nuclear 
shadowing and the cross section for the diffractive dissociation of the projectile j2oiJ23l|. which 



was established long time ago by Gribov p^. Its result can be derived using reggeon calculus 



24| and the Abramovsky-Gribov-Kancheli (AGK) cutting rules 25|] and is a manifestation 
of the unitarity. This formalism can be used to calculate directly cross sections of photon- 
nucleus scattering for the interaction with two nucleons in terms of the diffractive photon- 
nucleon cross section. In this formalism, the total photon-nucleus cross section is expressed 
as a series containing the contribution from multiple scatterings (1, 2, . . .): 

(1) I (2) I (3) I /n\ 

with the first term being the one that arises from independent scattering of the photon off 
A nucleons: 

a% = Aa^*p (4) 
and the first correction to the non-additivity of cross sections being 



a% = -47rA{A - 1) J d% Tlih) j^"^ dM' 



dM^dt 



FiiKin) (5) 

t=0 



where is the mass of the diffractively produced system. Fa is the nucleus form factor 
which takes into account the coherence effects and the differential 7*p cross section for 
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diffractive dissociation of the virtual photon appearing in (|5]) is given by: 



dM^dt 

where Bd is the diffractive slope parameter and X]pF!^{Q'^ ^xp, fi) is the diffractive pro- 
ton structure function. Moreover, tmin = — m^x^, xp = and (3 = Q'^/{Q'^ + M^). 
The integration limits in are M^j^ = 4m^ = 0.08 GeV^, M^^^ = {x^ra^^Jx - 1) 
and Xjpmax = 0.1. A shortcoming of this approach is that the inclusion of the higher order 
rescatterings is model dependent. This resummation is specially important at small x, where 
multiple scattering is more likely to happen. In general it is assumed that the intermediate 
states in the rescatterings have the same structure and two resummation schemes are con- 
sidered: (a) the Schwimmer equation 26|], which sums all fan diagrams with triple pomeron 
interactions and which is valid for the scattering of a small projectile on a large target. It 
implies that the photon-nucleus cross section is given by: 

r) = a,.,ix, r)Aj d'^Y^lA^Wm^) 
and (b) the eikonal unitarized cross section, given by 

a ,^{x, r) = a^.p(x, r) A d b — , (8) 



where 



As shown in 



2(A-l)/(x,Q2) 



Air /"^max drr'^ 

' xFj(t^i„). (9) 



f(x,Q') = -X dM 
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a^*p{x,r) JmI^^ dM^dt 



t=o 



23[, the eikonal unitarization predicts a larger magnitude for the nuclear 
shadowing than the Schwimmer equation. For models which take into account the possibility 
of different intermediate states see, e.g., Ref. j^^]. Except for the choice of the resummation 
scheme, the predictions for a^*A obtained using ([7]) or (|8]) are parameter free once the 
diffractive cross section is provided. Models based on this non-perturbative Regge-Gribov 
framework are quite successful in describing existing data on inclusive and diffractive ep and 



eA scattering 23 



|. However, they lack solid theoretical foundations within QCD. It is 
important to emphasize that some authors j2l| use these models as initial conditions for 
DGLAP evolution. 

The comparison among the predictions of the different models for nuclear shadowing 
presented in Ref. [ll, including the models discussed above, shows that they coincide within 
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~ 15% in the region where experimental data exist (x > 10~^) but differ strongly for smaller 
values of x, with the difference being almost of a factor 2 at x = 10~^. Our goal in this 
paper is try to reduce the theoretical uncertainty present in these predictions. In particular, 
differently from previous studies, which consider different inputs in the calculations using 
the Glauber, Schwimmer and Eikonal approaches, we will consider a unique model for the 
projectile - nucleon interaction. We will calculate the dipole - nucleon cross section and 
the diffractive structure function using the dipole picture and the solution of the running 
coupling Balitsky-Kovchegov equation 29|, which is the basic equation of the Color Glass 
Condensate formalism. Recently, this approach was shown to describe quite well the ep 
HERA data for inclusive and diffractive observables (See, e.g. Refs. 30|-|33|). Following this 



procedure we are able to estimate the magnitude of the theoretical uncertainty associated 
to the way the multiple scatterings are considered, reducing the contribution associated to 
the choice of initial conditions used in the calculations. Moreover, we discuss the possibility 
of discriminating between these unitarization procedures in a future electron-ion collider. 

This paper is organized as follows. In Sec. [Ill we present a brief description of inclusive 
and diffrative 7 - nucleon processes in the color dipole picture with particular emphasis 
in the dipole - proton cross section given by the Color Glass Condensate formalism. In 
Section IIIII we present the predictions of the three unitarization schemes discussed above 
using as input the CGC results for the dipole - proton interaction and compare them with 
the existing experimental data. Moreover, we present a comparison between the predictions 
for the kinematical region which will be probed in a future electron - ion collider. Finally, 
in Section [IV] we summarize our results and present our conclusions. 



II. INCLUSIVE AND DIFFRACTIVE 'jp PROCESSES IN THE COLOR DIPOLE 
PICTURE 

The photon-hadron interaction at high energy (small x) is usually described in the infinite 
momentum frame of the hadron in terms of the scattering of the photon off a sea quark, which 
is typically emitted by the small-x gluons in the proton. However, as already mentioned in 
the introduction, in order to describe inclusive and diffractive interactions and disentangle 
the small- a; dynamics of the hadron wavefunction, it is more adequate to consider the photon- 
hadron scattering in the dipole frame, in which most of the energy is carried by the hadron. 



6 



while the photon has just enough energy to dissociate into a quark-antiquark pair before the 
scattering. In this representation the probing projectile fluctuates into a quark-antiquark 
pair (a dipole) with transverse separation r long before the interaction, which then scatters 



off the target 



34| . The main motivation to use this color dipole approach is that it gives a 



simple unifled picture of inclusive and diffractive processes. In particular, in this approach 
the proton structure function is given in terms of the dipole - proton cross section, adp{x, r), 
as follows: 



d r dz\ilj{r,z)\ adp{x,r) 



(10) 



where 



is the probability of the photon to split into a qq pair of size r. Moreover, 



the total diffractive cross sections take the following form (See e.g. Ref. 35| ) 



1 da' 



V 



t=o 



Bd dt 



i=0 



where 



dt 



t=0 



^ j dh da\^T,L{a,r)\^<J%[ 



(11) 



(12) 



It is assumed that the dependence on the momentum transfer, t, factorizes and is given 
by an exponential with diffractive slope Br,. The diffractive processes can be analysed in 
more detail by studying the behaviour of the diffractive structure function F^^^\Q'^, (3^X]p). 



Following Ref. 



35| we assume that the diffractive structure function is given by 



(13) 



where the qqg contribution with longitudinal polarization is not present because it has 
no leading logarithm in Q^.^_The different contributions can be calculated and for the qq 



contributions they read [36|, l37| 



xjpF^,JQM3,xjp) 



32tt^13Bd 



X]ej2 f^^^ daa^{l-af%, 

r Jan 



(14) 



xpF'^^JQ'\l3,X]p) 



1287r^l3BD 



/•1/2 r 

Y,e}2 / daa{l - a) le'^la'^ + {1 - af]<l>i + m}<!>o} 



(15) 



where the lower limit of the integral over a is given by ao = | |^1 — ^ 1 — -^^ ) , the sum is 
performed over the quark flavors and 38 1 

2 



rdrKo^i{er)adp{x]p, r)\i{kr] 



(16) 



The qqg contribution, within the dipole picture at leading InQ^ accuracy, is given by |35l-l37| 

2 / \ 2 



z 



(17) 



ndn adp{u/kt,X]p)K2 



-\ 2 



'0 * \\i 1 — z 

As pointed in Ref. 39|], at small /3 and low Q^, the leading ln(l / /3) terms should be resummed 
and the above expression should be modified. However, as a description with the same 
quality using the Eq. ( |T7I) is possible by adjusting the coupling [39|, in what follows we will 
use this expression for our phenomenological studies. We use the standard notation for the 
variables xjp = (M^ + Q^)/{W^ + Q^) and x = Q^/{W^ + Q^) = /Sxp, where W the total 
energy of the 7*p system. 

The main input for the calculations of inclusive and diffractive observables in the dipole 
picture is adp{x,r) which is determined by the QCD dynamics at small x. In the eikonal 
approximation, it is given by: 



(Tdp(x,r) = (fbM{x,r,b) 



(18) 



where J\f{x,r,b) is the forward scattering amplitude for a dipole with size r = \ r \ and 



impact parameter b which can be related to the expectation value of a Wilson loop 12|. It 
encodes all the information about the hadronic scattering, and thus about the non-linear 
and quantum effects in the hadron wave function. In general, it is assumed that the impact 
parameter dependence of can be factorized as A/'(x, r, b) = A/'(x, r)S'(6), where S{b) is 
the profile function in impact parameter space, which implies adp{x,r) = (roAf{x,r). The 
brward scattering amplitude A/'(x, r) can be obtained by solving the BK evolution equation 



30| or considering phenomenological QCD inspired models to describe the interaction of 



the dipole with the target. BK equation is the simplest nonlinear evolution equation for the 
dipole-hadron scattering amplitude, being actually a mean field version of the first equation 



of the B- JIMWLK hierarchy 



40| . In its linear version, it corresponds to the Balitsky-Fadin- 



Kuraev-Lipatov (BFKL) equation [41|. The solution of the LO BK equation implies that 



the saturation scale grows much faster with increasing energy {Ql 



X 



with A ^ 0.5) 



than that extracted from phenomenology (A ~ 0.2 — 0.3). 

In the last years the next-to-leading order corrections to the BK equation were calculated 



42 



- |44| through the ressumation of agNf contributions to all orders, where Nf is the number 
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FIG. 1: Comparison between the predictions of the distinct models and the E665 experimental 
data at small x. 

of flavors. Thanks to these works it is now possible to estimate the soft gluon emission and 
running coupling corrections to the evolution kernel. The authors have found out that 
the dominant contributions come from the running coupling corrections, which allow us to 
determine the scale of the running coupling in the kernel. The solution of the improved BK 



equation was studied in detail in Ref. 43|. The running of the coupling reduces the speed 



of the evolution to values compatible with experimental data, with the geometric scaling 
regime being reached only at ultra-high energies. In 3o| a global analysis of the small x data 
for the proton structure function using the improved BK equation was performed (See also 



Ref. 



45|). In contrast to the BK equation at leading logarithmic asln{l/x) approximation, 



which fails to describe the HERA data, the inclusion of runnin g co upling effects in the 
evolution renders the BK equation compatible with them (See also 31n33|). In what follows 
we consider the BK predictions for ^/{x, r) (from now on called rcBK) obtained using the 
GBW \m initial condition. 
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III. NUMERICAL RESULTS AND DISCUSSION 



In what follows we shall consider two different nuclei, Ca and Pb, and use the deuteron 
(D) as a reference to calculate the experimentally measured ratios Rca/D = (2/40) ^2^"/^^ 
and Rpb/D = (2/208)^2^^/^2^. assume that the diffractive slope parameter is Bd = 6.7 
GeV~^ and that the nucleus form factor is given by: 



In Fig. [T]we compare the predictions of the Glauber (solid line), Schwimmer (dot-dashed 
line), Eikonal (dashed line) and double scattering (dot- dot- dashed line) models for the ratios 
with the E665 experimental data at small x Although joined with lines, our results are 
computed at the same (x) and (Q"^) as the experimental data. Our results demonstrate 
that if we compute the nuclear structure function up to two scatterings, which implies that 
(T^A = o'^M + CT^'A; we are not able to describe the experimental data. Furthermore, since 

(2) 

the magnitude of the first correction, a^.^^, is very large, then there is no hope to estimate 
the nuclear structure function by just summing a few terms in the multiple scattering series. 
Therefore, a full resummation of the multiple scatterings is necessary, which makes the 
predictions model dependent. The agreement with the current experimental data at small 
X of the Glauber, Schwimmer and Eikonal models is quite reasonable taking into account 
that no parameters have been fitted to reproduce the data. This implies that the current 
data are not able to discriminate between the unitarization schemes. 

Having in mind that a future electron - ion collider is expected to be able to analyse 
the kinematical region of small x (x ^ 10"^) and > 1 GeV^, we now compute the ratios 
Rca/D and Rpb/D SlS a. function of x for two different values of (= 1 and 10 GeV^). In Fig. 
[2] we present our predictions for = 1 GeV^. It is important to emphasize that in electron 
scattering the range of x- values attainable is kinematically restricted to x > Q^/s, where s 
is the squared center-of-mass energy, which implies that at = 1 GeV^ the smaller values 
of X in the perturbative region will be probed. At large x (~ 10"^) the predictions almost 




(19) 



where the thickness function is given in terms of the nuclear density pA as: 




with the normalization fixed by / d^b TA{b) = 1. 
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FIG. 2: Nuclear ratios Rp^^j^, (left panel) and Rca/D (right panel) as a function of x at = 1 
GeV2. 

coincide. However, at small x, the predictions based on the Schwimmer equation or on the 
eikonal unitarized cross section give a stronger shadowing than those based on Glauber-like 
rescatterings. In particular, at x ~ 10"^, the difference between Glauber and Schwimmer is 
almost 10 % in the ratio R{Ca/D) increasing to ^ 20 % in R{Pb/D). At this x value, the 
difference between Schwimmer and Eikonal is ~ 5 % and 12 % for the ratios R{Ca/D) and 
R(Pb/D), respectively. At smaller values of x, the difference between the three predictions 
increases, being larger than 20%. Consequently, a measurement of at A = Pb at small x 
with ^ 10% precision would be a sensitive test to discriminate between the different models. 

In Fig. [3] we present our predictions for the ratios Rca/D and Rpb/n as a function of x at 

= 10 GeV^. The behavior is similar to the one observed in the Fig. |2l The main point is 
that the differences between the predictions is not reduced significantly and this makes the 
discrimination between them possible also at this value of Q^. 

A final comment is in order. The results shown in Figs. [2] and |3] demonstrate that there 
is a large uncertainty associated to the choice of unitarization scheme used to treat the 
multiple scatterings and that, in principle, an experimental analysis of the nuclear ratios 
can be useful to discriminate between these approaches. 

Another uncertainty present in the study of the nuclear effects is related to the transition 
between the linear and nonlinear regimes of the QCD dynamics. We do not know precisely in 
which kinematical region the predictions obtained using the linear DGLAP evolution cease to 
be valid. In Fig. |4]we present a comparison of our predictions with those obtained using the 
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FIG. 3: Nuclear ratios Rpb/D (left panel) and Rca/D (right panel) as a function of x at = 10 
GeV2. 

EPS09 [l^ parametrization of the nuclear parton distribution functions, which is based on a 
global fit of the current nuclear data using the DGLAP dynamics. As it can be seen, due to 
the large theoretical uncertainty in the DGLAP prediction in the small-x region, represented 
by the shaded band in the figure, it is not possible to draw any firm conclusion about 
which is the correct framework to describe this observable in future eA colliders. This same 
conclusion was already obtained in [l^ in a somewhat different approach. Consequently, 
the study of other observables, such as the nuclear diffractive structure function 4a, and 

nn 

nuclear vector meson production pJl, |48| , should also be considered in order to discriminate 
between the linear and nonlinear regimes. To summarize: in order to learn more about the 
unitarization schemes using the nuclear ratios we must disentangle the nonlinear and linear 
regimes of the QCD dynamics. Our estimates show that due to the large freedom present in 
the DGLAP analysis they predict similar magnitudes for the nuclear ratios, which implies 
that a combined analysis of several observables is necessary. 



IV. CONCLUSION 



The behaviour of the nuclear wave function at high energies provides fundamental infor- 
mation for the determination of the initial conditions in heavy ion collisions and particle 
production in collisions involving nuclei. One of the main uncertainties is associated to the 
magnitude of the nuclear shadowing, which comes mainly from the way in which the multiple 
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FIG. 4: Predictions of the different models discussed in the text. The dash-dash-dot hne repre- 
sents the central value of the prediction obtained with the EPS09 parametrization of the nuclear 
parton distribution functions. The shaded band represents the theoretical error coming from the 
uncertainties in the EPS09 parametrization. 

scattering problem is treated and from the modelling of the projectile - nucleon interaction. 
Since a future EIC will probe the shadowing region while keeping sufficiently large Q^, new 
studies which determine the main sources of uncertainties in the predictions are necessary. 
In this work we compare three frequently used approaches to estimate the nuclear shadowing 
in nuclear DIS. As in these approaches the nuclear cross section is completely determined 
once the interaction of the projectile with the nucleon is specified, we considered a single 
model (rcBK) as input of our calculations in order to quantify the theoretical uncertainty 
which comes from the choice of the unitarization model. In particular, we calculate the 
nuclear ratio between structure functions considering the Glauber, Schwimmer and Eikonal 
approaches down to very low-x utilizing the rcBK results both for inclusive and diffractive 
cross sections in 7*p scattering. Our results demonstrate that the current experimental data 
at small x are described successfully by the three approaches. However, the difference be- 
tween their predictions becomes large in the kinematical region which will be probed in the 
future electron - ion colliders. 
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